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Abstract 

o 

This note collects formulas relative to the energy dependence of the cross sections for open charm 
and J/ip production for iViV-collisions at threshold. It is a basis for the best input to MonteCarlo 
calculations, for associative charm particle production in nucleon-nucleon, nucleon- nucleus, ion-ion 
and proton-antiproton collisions. 
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FIG. 1: Feynman diagram for p + p — ► A+ + D + p. 

I. THE REACTION p+p^A++F + p 

The threshold behavior of the total cross section for the reaction p + p — > A+ + D° + p 
can be described by: 

a(pp -> A+^Op) ~ 0.2(Q/0.1 GeV) 2 nb. (1) 

with 

Q = ^i- (m + M + j/), (2) 

m is the nucleon mass, M{p) is the mass of A C (D). This estimation is based on D— exchange 
mechanism, Fig. [TJ and on the comparison with the processes p + p — > A + + p 
(if— exchange), assuming SU(A) symmetry [lj]. The threshold Q 2 -behavior ( Q is the en- 
ergy excess over threshold) ,is the phase-space for a three-body reaction 

In Fig. [2] we report the threshold behavior of the cross section for the reaction p + p —>■ 
A+ + D° + p as a function of Q, (top) and of the proton kinetic energy (bottom). 

No experimental data exist in the threshold region. The lowest energy where the open 
charm cross section has been measured, E p =70 GeV, yfs = 11.46 GeV, corresponds to Q = 6 
GeV. From Eq. (CQ) we find, for this energy: a{pp — > A+D°p) ~ 0.7 fib, to be compared to 
the experimental value a(pp — > charm) = 1.6iJ;y (siat) ±0.3 (syst) fib [3]. The calculation 
rigorously hold only for Q < \i. 
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FIG. 2: Cross section for the reaction p + p — > A+ + D° + p as a function of the energy excess over 
threshold, Q, (top) and of the proton kinetic energy (bottom) 

II. TOTAL CROSS SECTION FOR OPEN CHARM PRODUCTION 



The total cross section for pp-interaction near threshold can be approximated as: 



a tot (pp) ~ 40 mb at yj(s) > 10 GeV 

and it flattens at larger energies. 

The semi- inclusive cross section for charm particle production has been calculated in jlj] , 
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in analogy with strange particle production, and it has been found 

a(pp -> XA C ) ~ 40 fib 

depending on the choice of the cut off parameter for the form factor, which has to be 
introduced in the vertex of the diagram. Let us note that this number is larger than the 
direct experimental value Using these cross section, one can find: 

R D = aiPP ^ AcX) = (0.25 -hi) -HT 3 

Following [5j, in analogy with strange particle production in 7iV-collisions, in principle 
one should take into account another contribution, D*-exchange, which may increase the 
cross section of a factor four. Therefore: 

R D+D * ~ (l-4)-l(r 3 

which is in agreement with the experimental value of 40-r- 200 /ib at ^=62 GeV Q]. 



A. Isotopic relations 

The models allow to calculate the cross section for all channels of exclusive charm 
particle production, not only for pp- but also for np-collisions. The starting point of such 
calculations is the isotopic structure of the matrix element M. for iViV — > Y c + D + N, which 
is model independent. 

We give here the necessary relations based on the isotopic invariance of the strong inter- 
(£„) = 1, I (A c ) = 0, 1(D) = 1/2: 



action, with / 
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where Aj x i 2 are the isotopic amplitudes, corresponding to the total isospin l\ for the initial 
nucleons and the total isospin I2 for the produced DiV-system. In fact we have twice more 
of these reactions, by changing isospin p <-> n (for example pp — > T^ + D^p <->• nn — > 

One can see, from Eq. (jHJ), that seven different matrix elements for different iV + N — > 
S c + D + A^-reactions, are characterized by three complex isotopic amplitudes: An, A\o and 
Aqi. Relations flSJ) hold for any model in any kinematical condition. 

Model independent relations among the seven differential cross sections hold: 



01 + 02= 2\A n \ 2 +4|Ai | 2 
^3= 2\A U \ 2 , 



a 4 + a 5 = 2\A U \ 2 +4|A i 

loi 



(4) 



a- 6 + (7 7 = 2|A 10 | 2 +2\A 12 



The isotopic averaged cross section for the processes iV + iV — > £ c + D + A^, typically 
used in the transport codes, is defined as: 

1 



Therefore: 



a = 7 • 2 • E °i 

^ i=l-7 



o = 3(|A n | 2 + |A 10 | 2 + |A 121 



oil 



and does not contain any interference term between the different isotopic amplitudes. The 
isotopic amplitudes can be expressed in terms of the cross sections as: 

l^ll| 2 = |*3, 

I An I 2 = |(o- 4 + 0-5 - 0-3) (5) 
l^iol 2 = + °>) - 7(04 + 0-5 - 0-3) 

The situation essentially changes near the reaction threshold, where all final particles 
are produced in 5-state. The selection rules with respect to the Pauli principle, the P- 
invariance and the conservation of isospin and of the total angular momentum, allow to 
find a more simple spin structure for the matrix element for pp and np-collisions in Q. A 
direct consequence of the threshold spin structure is that the interferences An (g) A^ L and 
Axq® Aq X , which can be present in the general case, vanish in the threshold region, after 
summing over the polarizations of the colliding nucleons. Therefore at threshold we can 
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derive the following relations, additional to (jl]), which hold for any model: 

^{np -> E++ZT™) = ^-[np -> E^p), 
cto; au; 

-^-{np — »■ T.^D^p) = —{np — > £+.D ?7,). 



B. Scheme for estimation of cross sections 

Taking the following assumptions: 

1. D-exchange model: the amplitudes are proportional to the amplitude for elastic DN- 
scattering, which is a combination of amplitudes for the two possible states of isotopic 
spin, A\ and Aq. 

2. SU(4) symmetry: A\ 3> Aq (as for KN scattering), so we neglect the amplitudes 
corresponding to I(DN) = 0. 

3. Threshold regime: all particles are in S final state; no interference between singlet and 
triplet initial states. 

We find: o\ = <7 2 = -a 3 = |A n | 2 , cr 4 = a 5 = \Au\ 2 + 2|A i| 2 = 3a\, and er 6 = a 7 = 
|A)i| 2 = 01, cr 9 = a w = l/4a 8 . 

This derives from the fact that An and Aoi are both proportional to Ai, therefore the 
spin structure of theses amplitudes may be different, but their moduli squared are the same. 



<7i can be calculated from the proportionality with a 8 = a(pp — > A^D°p): a x = | c ag. 

9pK c D 

If one assume S77(4)-symmetry, in the considered D-exchange model, the ratio of the 

9 2 t, d q\k 
coupling constants | c can be related to -| — . An estimation of this ratio from the 

9pK c D 9 P AK 

literature (j^j and refs herein) gives ^ ScD ~ 0.1. 

9 p k c d 

The overall D°-production, due to the discussed reactions, results in about 2cr 8 . 

C. Results 

The calculation [1] takes into account the general threshold symmetry properties of the 
strong interaction, such as the the conservation of parity, total angular momentum, and the 
Pauli principle. 



Current models, which reproduce high energy data, can not always be safely extrapolated 
to threshold. For example, in J] the cross section is taken as: 

a {PP ^ A+P°p) =a x (l- Z) a Z~P (6) 



with a x = 0.496, Z = yjs° x /^, ^ s° x = 5.069 GeV, a = 4.96 and (3 = 1.36 for D° 
production, and slightly different paramaters for D "-production: ax = 0.363, \J~Sx~ = 5.073 
GeV, a = 4.94 and (5 = 1.44. 

In Fig. [3] a comparison between the calculations [1] (thick lines) and 4] (thin lines) is 
shown for D° production (solid lines) and for D~ -pro duct ion (dashed lines). 

One can see that the values can differ by an order of magnitude. This comes from the 
fact that the correct Q-behavior at threshold, from phase space considerations, should be 



quadratic, whereas the best fit parameters of Ref. [4] give an exponent ~ 5. The difference 
is also due to the fact that the isotopic relations among D° and /}~-cross sections are not 



respected in the parametrization [4]. This induces very large effect in particular in the 
threshold region. 

Note that the calculation which is based on the Quark-Gluon String model and the 
Regge phenomenology, is in agreement with [l| , for a specific choice of form factors. 

D. Possible future studies 

One can foresee an extended program for the study of fundamental processes of charm 
hadropro duct ion, which will be useful for the experiments at the future international facility 
at GSI: 

1. Analysis of vector charm Z)*-associative production in AW-collisions N + N — > Y c + 

D* + N; 

2. Study of associative production of hperonic and nucleonic resonances: N + N —>■ 
Y c * + D + N and N + N -> Y c + D + N*; 

3. Study of associative charm production in AN, A A and BB collisions, where BB is 
any baryonic or hyperonic resonance, which are important for high energy ion-ion 
collision; 

4. Study of the transition regime for threshold to Regge description; 
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FIG. 3: Comparison between the calculations 



1] (thick lines) and 



(thin lines) for D° production 



(solid lines) and for D -production (dashed lines). 



5. Consideration of diffractive production of charm particles in pp and pp-collisions. 



III. CROSS SECTION FOR p + p -> p + p + J/ip AT THRESHOLD 

A general formalism for the calculation of the cross section and the polarization observ- 
ables for the reaction p + p^p + p+ J/ip at threshold has been developped in [9|. From 
a model independent parametrization of the matrix element, valid for all types of vector 
meson production in NN collisions, it is possible to derive expressions for the polarization 
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phenomena and for the ratio of J /if) production in np and pp collisions, in terms of the 
threshold amplitudes. In framework of a model based on t-channel meson exchange and 
comparing with the existing data on production, it is possible to predict the threshold 
behavior of the cross section. 

The ratio R of cross sections for and J /if)- production in framework of the same ap- 
proach, namely for 7r-exchange in N + N — > N + N + V° and p-exchange for the subprocess 
n + N — > iV + V°, with V° = <f> or J/if>, one can find: 



a(pp ^ ppJ/ip) g 2 (J/ip -> np) ( t^-ml\ 2 f t 4 



m 



a(pp -> ppfy g 2 (<P^np) \t.J/4>-mlJ \tj/^-m 



V ps (J/i>) F{<j>) 



v ps (<p) HJ/^P) 

where g(V — > np) is the coupling constant for the decay V — > np, t v = —mmy is the 
threshold value of the momentum transfer squared, F(t) is a phenomenological form factor 
for the vertex n*p*V°, with virtual n and p. 

Taking the volume of three-particle phase space as 2j, one can find: 



V ps {4>) M \2m + M J/lP ) 

and for the particle flux: 

T(<fy _ k{(j>) ( 2m + M 4> 



3/2 



T(J/if)) k{Jli>)\2m + M m )' 

where k is the momentum of the initial nucleon in the CM system: k = \J~^~ m ^ s = 
(2m + my) 2 . 

Using the existing experimental data about the decays J /if) — > n + p and <f> — > n + p , one 

can find g 2 (J/ip — > np)/g 2 (<p — > np) ~ 10~ 4 , 

( + —m 2 \ 2 ( tA, — m 2 \ 2 

Taking the terms: ^ =0.111 GeV 2 and — p — =0.196, the ratio of 

ytj/^-mlJ VJ/i>- m2 P ) 

phase space: ps ^ — %- ^.^} . = 0.216, one finds for the ratio of cross sections, at the same 
V ps {<f>) F{J/ifj) 

Q: 

ajjpp ^ppj/ip) = 4 7 _ 1Q _ 7 
a(pp — > pp0) 
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The relevant available experimental data are: a(pp — > pp<fi) — 300 nb at pi = 3.67 GeV 



10|, and a exp (pp ppJ/ib) 



One can find that 



0.3 ± 0.09 nb for ^ = 6.7 GeV (Q=1.725 GeV) [8|. 



a(pp — > ~ 206 



0.1 GeV 



nb 



and 



a{pp — » ppj/ip) ~ 9.7 ■ 10 



-5 



F{tj,^)/Fm) 



nb . 



(7) 



.0.1 GeV 

This value is too small, when compared with the existing experimental value for the 
lowest <Js = 6.7 GeV, namely a exp {pp — > ppJ/ip) = 0.3 ± 0.09 nb. 

Note that the p-exchange model for a(irN — > J lib) gives a cross section one order of 

nnn 

magnitude smaller in comparison with other possible theoretical approaches pjj, 1121, H3J . It 
is one possibility to explain the value of a exp {pp ppJ /ip). Another possibility is to take 



F{tj^)/F{^) 



10, which can be plausible, because the J/tb — cc-system must have a 



smaller size in comparison with 



ss. This can be realized by the following form factor: 
1 



Fy{t) 



t 



with Ay ~ my. 

The cross section, based on Eq. ([7]) normalized the the experimental point at ^/s = 

r -I 2 

6.7 GeV, i.e., taking the ratio F(tj/^)/ Fit^) = 10, is plotted in Fig. [7J together with 

the experimental data from the compilation 14], where different symbols differentiate Jib 

production in pp or extrapolated from pA collisions. Note, in conclusions, that in the 

framework of the considered model, one can find: 

a{np — > npj/ip) 
a{pp — > ppJ /ip) 

which would require a correction of the experimental data on pA reaction, where equal np 
and pp cross sections are usually assumed for the extraction of the elementary cross section. 
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FIG. 4: Cross section for J /tp production in pp collisions. Data are 
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